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NYO-9600 
Magnetohydrodynamic  Nozzle  Flow  with  Three  Transitions 

Abstract 

The  flow  of  an  infinitely-conducting  inviscid  gas  in  the 
presence  of  a  magnetic  field  everywhere  parallel  to  the 
velocity  is  considered.   Such  a  flow  starting  from  rest  and 
expanding  to  supersonic  velocities  has  three  transitions, 
across  each  of  which  the  flow  differential  equations  change 
type:   from  elliptic  to  hyperbolic  or  vice  versa.   The  flow 

velocities  at  the  three  transitions  are,  respectively,  the 

P    PI  /p 
sound  speed  c,  the  Alfven  speed   a,   and  the  speed  ca/(c  +  a  )  ^ 

A  class  of  those  flows  is  approximately  constructed  by  solving 
the  flow  differential  equations,  with  the  flow  velocity  etc., 
prescribed  on  an  axis  of  symmetry.   The  solution  is  in  the 
form  of  a  formal  power  series  in  a  neighborhood  of  the  axis. 
The  resulting  streamlines  give  possible  symmetric  nozzle 
contours.   Some  mathematical  investigations  are  made  to  de- 
termine the  acceptable  form  of  the  data  to  be  given  on  the 
axis.   Numerical  examples  giving  nozzle  shapes  are  included 
and  compared  with  similar  nozzles  for  non-conducting  gases. 
The  procedure  described  is  an  extension  of  that  of  Friedrichs 
for  ordinary  gas  nozzles. 
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1,   INTRODUCTION 

This  paper  gives  a  formal  power  series  solution  (near  the 
axis  of  symmetry)  for  the  plane  nozzle  flow  of  an  infinitely- 
conducting  inviscid  gas  with  a  magnetic  field  everywhere  parallel 
to  the  flow  velocity.   The  gas  expands  from  near-reservoir  condi- 
tions to  supersonic  velocities  in  the  nozzle. 

When  an  ordinary  non-conducting  gas  expands  from  rest  to 
supersonic  velocities  in  a  de  Laval  nozzle,  it  is  well-known 
that  the  gas  crosses  a  sonic  transition  (near  the  throat  of  the 
nozzle) J  upstream  of  which  the  flow  is  subsonic  and  governed  by 
an  elliptic  differential  equation,  and  downstream  of  which  the 
flow  is  supersonic  and  governed  by  a  hyperbolic  differential 
equation.   Moreover,  we  know  that  such  a  flow  is  stable,  whereas 
the  reverse  flow  (i.e,  supersonic  flow  slowing  down  to  subsonic 
flow  through  continuous  compression)  is  unstable,  and  thus  occurs 
physically  only  when  accompanied  by  shocks. 

In  the  case  of  an  infinitely-conducting  gas  with  a  magnetic 
field  parallel  to  the  velocity  everywhere,  it  is  by  now  well- 
known  (See  e.g.  Grad  )  that  if  such  an  expanding  flow  is  poss- 
ible, the  gas  must  cross  three  transitions.   The  differential 
equation  governing  the  flow  is  elliptic  at  low  velocities, 
changes  to  hyperbolic  across  the  first  transition,  back  to 
elliptic  across  the  second,  and  again  to  hyperbolic  across  the 
third.   For  def initeness ,  we  shall  assume  that  the  field 
strength  is  such  that  the  gas  reaches  Alfven  speed  "a"  before 
it  reaches  sound  speed  "c"  (although  the  procedure  is  general 
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so  that  assuming  the  converse  is  just  as  good).   The  first 

2   2  1  /2 
transition  occurs  at  a  gas  speed  equal  to  ca/(c  +a   )  , 

the  second  at  "a"\    and  the  third  at  "c".   Thus  far^  we  have 

no  definite  knowledge  that  such  a  fully  "transonic"  flow 

exists.   Intuition  is  insufficient*  for  example,  Grad  has 

shown  that  the  addition  of  a  small  transverse  field  component 

will  prevent  such  flows  from  crossing  the  sonic  transition, 

a  fact  which  is  intuitively  non-obvious.   The  present  paper 

does  not  prove  the  existence  or  inexistence  of  such  fully 

"transonic"  flows,  nor  does  it  answer  the  question  whether 

such  flows  are  stable.   It  merely  gives  an  approximate 

description  of  such  flows  in  nozzles,  their  existence  being 

assumed.   Nevertheless,  the  well-behaved  properties  of  the 

solution  tend  to  suggest  the  probable  existence  of  these 

flows,  rather  than  their  inexistence. 

2 
The  procedure  used  is  that  of  Friedrichs   for  ordinary 

gas  nozzles.   The  potential  and  stream  functions  are  used 
as  independent  variables,  and  the  space  coordinates  and 
hodograph  variables  as  dependent  variables.   Data  are  given 
analytically  on  a  line  (say,  the  x-axis)  to  give  solutions 
symmetric  with  respect  to  it.   This  is  a  Cauchy  problem,  and 
the  solution  is  sought  in  the  form  of  a  power  series,  its 
existence  being  guaranteed  by  the  Cauchy-Kowalewski  Theorem, 
Any  streamline  in  the  physical  plane  is  then  a  nozz]e  con= 
tour.   Since  the  region  of  convergence  of  the  series  solu- 
tion is  unknown,  this  solution  is  only  a  formal  solution 


for  small  distances  from  the  axis  of  symmetry,  with  the 

first  term  equal  to  the  one-dimensional  or  "hydraulic" 

3 
solution.   By  a  similar  procedure,  Frankl  went  further  to 

investigate  the  region  of  convergence,  thus  yielding  an 
existence  proof  in  addition  to  an  approximate  description 
for  the  flow. 

Our  procedure  is  essentially  the  same  as  that  of 
Friedrichs.   We  assign  a  velocity  distribution  on  the  axis 
and  obtain  the  nozzle  contour  as  well  as  all  flow  variables 
in  the  neighborhood  of  the  axiSo   However,  as  will  be  seen, 
our  data  cannot  be  given  as  arbitrary  analytic  functions, 
because  the  axis  becomes  characteristic  at  one  of  the  tran- 
sitions» 

Section  2  gives  the  basic  differential  equations  of 
the  floWj  and  the  transformation  of  variables  that  are  used. 
Section  3  gives  the  formal  power  series  solution,  in  terms 
of  given  Cauchy  data  on  the  axis.   Section  ^■   is  devoted  to 
an  investigation  of  the  singularities  and  the  necessary  re- 
strictions on  the  Cauchy  data.   Finally ^  we  give  some  numer- 
ical  examples  in  Section  3»    and  discuss  the  results  and  con- 
clusions,, 

£,.  THE  BASIC  DIFFERENTIAL  EQUATIONS 

The  basic  differential  equations  for  the  steady  flow 
of  an  infinitely  conducting,  mviscid,  compressible  fluid 
(i.e.,  the  steady  state  Lundquist  equations)  follow: 

_  4  - 


Continuity    <\iv    (pq)   =   0  (1) 


Momentum     -^^  +  (q- V)a  +  -^  x  curl  B   =   0         (2) 


Energy       {q.V)S   =    0  (3) 

Equation  of  State      p  =   C(S)  p^  (4) 

Maxwell's  Equations    curl  (  B  x  q  )   =   0             (5) 

div  B    =  0              '    (6) 

Here  p   is  the  density,   q  the  velocity,   p  the  pressure, 
B  the  magnetic  field,  S  the  entr-opy,  and  y   t^©  ratio  of 
specific  heats.   This  is  a  system  of  10  equations  in  9 
unknovms .   Had  we  considered  the  time-dependent  Lundquist 
equations,  (6)  would  merely  play  the  role  of  an  admissi- 
bility condition  on  the  initial  data,  (as  is  well  known  in 
Maxwell's  equations),  and  (1)  to  (3)  form  a  well-defined 
hyperbolic  system.   In  steady  flow,  (6)  is  needed  as  an 
additional  equation.   We  consider  the  special  class  of 
such  flows  (see  e.g.  Grad  ,  Kogan-^,  Resler-Sears")  where 
the  field  and  velocity  are  everywhere  parallel' .   As  is 
done  in  all  these  references,  we  let  B  =  a  p  q.    Then  (5) 


would  be  trivially  satisfied.   (1)  and  (6)  require  that  a 
be  chosen  constant  along  each  streamline.   Moreover,  along 
each  streamline,  (2)  integrates  into  an  ordinary  Bernoulli 
equation 

^  +  -3^   =     E  (7) 

P      ^^ 

where  E,  the  total  enthalpy,  is  constant  along  each  stream- 
line. 

In  this  paper,  we  make  the  further  restriction  that 
each  of  the  three  quantities  S,  a,  and  E  be  constant  not 
only  along  every  streamline,  but  also  throughout  the  flow. 
This  would  be  the  case,  of  course,  if  at  infinity,  a  uni- 
form flow  with  a  parallel  field  exists.   Otherwise,  this 
is  an  additional  restriction. 

Rewriting  the  momentum  equation,  we  have 

V  E  -  q  X  curl  q  +  (a/M.)q  x  curl  apq  =  0 

—  p    — 

or      q  X  curl  (l  -  a  p/M.)q   =   0 

This  implies  that  curl(l-<i  p/M-)q   must  vanish,  if  it 
vanishes  at  infinity.   The  equations 

div  (pq)       =   0  (c) 

curl   (  1  -  a  p/m-  )  q   =   0 

are  indeed  the  exact  counterparts  of  the  continuity  and 
irrotationality  equations  in  ordinary  fluid  dynamics.   The 
quantity  a  p/jj.  is  readily  identified  as   (B  /M-p)/q  =  lA^  > 
or  the  reciprocal  of  the  Alfven  number  squared.   The  same 
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equations  could  be  derived  if  we  had  started  off  with  the 

o 

hydrodynamic  analog  of  Imai  ,  and  observed  that  if  a   is 
constant,  a  Kelvin's  Theorem  would  hold;  then  (8)  is  deduced 
in  the  same  way  as  the  continuity  and  irrotationality  equa- 
tions are  deduced  in  gas  dynamics  from  Euler's  equations 
using  the  condition  of  isentropy, 

V/e  now  consider  plane  flows",  and  introduce  a  potential 
function  cp   and  a  stream  function  i/        in  the  usual  way. 
Letting  q  denote  the  magnitude  of  the  velocity,  and   9   its 
inclination  to  the  x-axis,  we  have 

1^   =    (  1  -  1/a2)  q  cos© 

II   =    (  1  -  1/a2)  q  sin©  (^j 

^        =   -  P  q  sin© 


6y 


=    p  q  cos  0 

V/e  shall  consider  flows  which  are  symmetric  with 
respect  to  the  x-axis.   In  what  follows,  the  subscript  (  ) 
denotes  quantities  on  the  x-axis  (as  functions  of  x  alone). 
The  subscript  (  )j   denotes  quantities  at  the  first  transition 
(M'^+A-  =  q^/(a^c^/a^+c^)  =  1),  the  subscript  (  ) -^j   denotes 
quantities  at  the  Alfven  (in  our  case,  the  second)  transition 
(A  =  q/a  =  1),  and  the  subscript  (  )ttt  denotes  quantities 
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at  the  sonic  (in  our  case,  the  last)  transition 
(M  =  q/c  =1).   We  shall  consider  9  and  t  ,  or  rather, 
convenient  functions  of  each,  as  independent  variables; 
this  permits  the  flow  region  being  considered,  namely  in  a 
plane  nozzle,  to  be  mapped  into  a  rectangle „   (Fig.  1). 
To  this  end,  we  introduce  ?,    =   £^9 )  ,  ^l  =  n  ( ,|, )   thus: 

cp   =     (  1  -  l/j\^{x)    )  q^(x)  dx 

Jo  °        °  (10) 

*   =   PJII  ^III  "^ 

It  is  clear  from  the  first  two  equations  of  (9)  that 
if  9  and  t  were  used  as  independent  variables  directly, 
the  entire  flow  region  in  the  x-y  plane  would  be  mapped  into 
the  left  half  of  the  9-^  -plane  covered  twice,  with  a  fold 
at  A  =  1.   Introducing  i    avoids  this,  and  makes  the  mapping 
one-to-one.   t]  is  introduced  merely  for  the  convenience  of 
having  both  independent  variables  as  distances.   Physically, 
I    is  the  distance  along  the  x-axis  to  an  equipotential 
line  9  =  constant,  and  "H  is  the  distance  along  the  sonic 
line  (M  =  1)  to  a  streamline  ^      =  constant.   With  r,  and  1 
as  independent  variables,  the  equations  for  q,  0,  x,  and  y 
are: 

_A_    f  (    1    -    l/An^j    q^l     _  do  (^^) 

6n    I  (  1  -  i/A  2)  q-  ]   -    -^M  ^     ' 

(     1    -    l/Ao2)    go)    ^  ^      ^  ^^2) 

(    1    -    1/A    2)    q     I    <^^  ^^ 
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^,  =   -  h  sin  Q  (13) 

O  'I 


^     =  h  cos  0  (1^0 

c'l 

where  h  is  defined  as   h  =  p   q-nr/  pq   i.e.,  an  elementary 
area  ratio. 

The  equ-ttionn  (11)  -  (1^)  are  obtained  in  the  usual 
way,  by,  e.g.,  inverting  (9)  and  cross-differentiating. 
They  serve  as  our  basic  set  of  differential  equations. 
It  will  be  seen  in  Section  ^■   that  the  coefficients  are 
analytic  functions  of  i    ,  'H  ,   q,  and  Q  ,  thus  tending  to 
justify  the  use  of  the  Cauchy-Kowalewski  procedure  near 
the  X-axis. 

3.   POWER  SERIES  EXPANSIONS 

We  seek  a  solution  of  (11)  -  (14)  in  the  neighborhood 
of  the  axis  of  symmetry  by  posing  Cauchy  data  for  q  (s  ), 
Q  Ai) ,    X  (|),  and  yAi)    on  the  axis  ^  =  0,   Evidently, 

^0=  °'  ^o  =  ^  '  ^o  "-  °» 

As  the  equations  (11)  -  (14)  are  alternately  elliptic 
and  hyperbolic,  the  Cauchy  problem  is  not  well-posed  unless 
q  (^)  is  analytic.   The  existence  of  a  solution  in  some 
neighborhood  of  the  axis  is  then  guaranteed  by  the  Cauchy- 
Kowalewski  theorem,  provided  the  axis  is  everywhere  non- 
characteristic.   We  shall  not  investigate  the  region  of 
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convergence;  otherwise,  we  can  indeed  attempt  to  prove 
the  existence  of  the  solution  to  (11)  -  (l^)  in  this 
manner.   In  our  case,  we  shall  content  ourselves  to 
obtain  a  formal  series  expansion,  the  first  term  of  which 
will  turn  out  to  give  the  "hydraulic  approximation" .   At 
first  sight,  we  might  be  concerned  with  the  stability  of 
the  flow,  since  small  (non-analytic)  variations  of  the 
data  on  the  axis  can  produce  large  variations  in  the  solu- 
tion.  This,  however,  does  not  at  all  imply  the  instability 
of  the  flow;  for  in  the  actual  flow,  the  walls  of  the  nozzle 
will  be  fixed,  and  the  small  flow  perturbations  do  not  cor- 
respond to  perturbations  in  the  Cauchy  data  on  the  axis. 
This  method  is  valid,  therefore,  in  yielding  the  contours 
of  the  flow  nozzle,  which  would  be  any  of  the  streamlines 

in  the  x-y  plane.   Thus  far,  the  method  described  is  identi- 

2 
cal  to  the  Friedrichs  procedure   for  ordinary  gas  nozzles. 

A  rather  severe  difficulty,  however,  arises  at  the 

2   2 
first  transition  M  +A   =1.   Here,  the  axis  of  symmetry  is  in 

the  characteristic  direction.   At  A  =  1 ,  similar  considera- 
tions arise.   For  the  time  being,  however,  we  shall  ignore 
those  questions  and  merely  obtain  the  formal  series,  leaving 
until  Section  4  a  more  detailed  investigation  to  see  how  the 
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initial  data  may,  and  must,  be  restricted  at  these  singu- 
larities. 

We  shall  naturally  give,  as  known  data,  the  three 
constants  S,  a  ,  and  E,  and  the  stagnation  density  'p'  , 
On  account  of  (7),  all  the  ordinary  isentropic  gas- 
dynamical  relations  (between  velocity,  Mach  number,  and 
the  thermodynamic  variables)  obtain.   In  particular,  the 
values  of  all  these  variables  at  the  three  transitions  can 
be  determined.   We  list  these  in  terms  of  the  given  data, 
and  shall  henceforth  treat  them  as  known  data  also. 

Sonic  transition: 


(15) 


%  =  1 


"^m  =  ^m    =    2(^^)  E 


p^   =   [2/(Y-l)]l/^Y-l)  - 


Alfven  transition: 


%       =      1 

p-jj       =     M/a^  (16) 

2 


=    (2/Y-l)     [(a2$/M)Y-l   _   i] 
I      =    (Y-1)    E    (u/a2^)Y-l 
q^      =  ^n     Cjj2      =      2E     [l    -    (M/a2^)Y-l  J 


M 


^n 
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First  transition: 


Mj^-  +   Aj^ 


Mj2  i3  given  by   [l  +  ( y-1)Mj^/2  ]  ^^^""^  =(a2p'/n)(l  -  M^^) 
Pj   =    P  [l  +  (y-1)m//2]  ^/^-Y  (17) 

Cj2  =    (y_i)  e  [l  +  (y-1)Mj2/2]  -1 

Let  now  q  (£),  hence  also  p  ,  h  ,  N   ,  etc.,  be  given 
analytically  in  2  .   Then  we  expand  in  r\      the  variables 
q,  0,  X,  and  y,  depending  on  whether  they  are  even  or  odd 
in  r\      (and  with  coefficients  functions  of  |   to  be  found), 
thus: 

q(^,Ti)   =   Qq^  ^  ""  ^2^^    ""  ^^^l^  +....) 
0(e,Tl)   =   0-^Ti   +  0^  n-  +   

,^       \             f           ^             2  ^     4  ^  (18) 

x(|,Ti)   =   ^       ^2  ^      '^ij-''^     


Now 


y(l,Ti)  =  y^Ti  +  y^  Ti^  +  

2 
p/p^    =   (cVc/)  1/Y-l    ^  1  .  M^2p^^2 


^n'^,     ^o   .  .  2    ..  ..^M  2  o  2v   ^ 

n  I  2 


2-(  2P^  +  ^2   -  (2-y)Mo^  ^2"^  "1   "  '" 

2  o 


h/h^    =   1  +  (m/-1)P2^2  -  l^^l-^"^)  -  ^P2^D-^2-Y)M^] 

(19) 


-  (1  -  M^^)^32^|  "^^  - 
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-1 


(1-1/Ao  )qo   ^  (£  ^"   f  p/po  -  Aq 
(1-1/A2  )q 


-1 


=  1  - 


-  3, 


(-^j.3^^-|^(^^ 


n 


o  o  o  «  O  O 


Substitution  of  (18)  and  (19)  into  (11)  -  (1^)  yields 
for  the  0-th  and  1st  powers  in     , 


0,  =  dh^/d^  =  h^' 
1      o        o 


(20) 


This  is  evidently  the  "hydraulic  approximation",  and 
it  is  in  no  way  different  from  ordinary  gas  dynamics,  as 
has  been  noted  by  Grad  .   The  next  higher  order  terms  give 


h  = 


JloJln"'  An^  -  1 


2   IMq2-,a/-  1/ 


^3  =  3  6r 


hn^hn"      (Mn2-l)(An^-l)]    ^   hoho'hn" 


\ 


M 


^+A  "".    1 
o       o      / 


^3  =  ''o 


(M^2_i)3^/  3   _   Q_^2/  ^ 


~i 


(21) 


^2  =  -\   ^'  /  2 


(20)  and  (21)  give  the  first  and  second  order  coeffici. 

ents  for  the  expansions  of  q,  0,  x,  and  y  (equation  (18)  ) 

near  the  axis  of  the  nozzle. 
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Evidently,  h  may  not  be  arbitrarily  assigned  near  the 
first  transition,  for  p^  will,  in  general,  become  infinite. 
This  is  not  surprising  since  the  axis  is  characteristic 
there  J  as  was  pointed  oui-  earlier.   To  see  how  this  initial 
data  must  be  restrictedj  we  attempted  to  continue  the  process 
to  obtain  the  higher  order  coefficients,  i.e,  3^.;  P^   etc., 
and  to  observe  from  these  some  sufficiency  requirement  for 
the  p's  to  remain  finite.   For  example,  from  the  expression 
of   gp  in  (21);  h  "  must  evidently  be  required  to  vanish  at 
the  transition.   This  procedure j  however,  led  to  some  rather 
inconclusive  results|  the  only  obvious  condition  thus  deduced 
was  that  all  higher  order  derivatives  of  h  .  from  the  second 
onwards,  should  vanish.   Clearly,  this  is  too  restrictive. 
It  results  m  a  linear  distribution  of  h  .  which  is  merely 
a  source  or  sink   Needless  to  say,  the  source  or  sink  flow 
cannot  cross  all  three  transitions,  since  M  =  1  is  a  limit- 
ing line  for  it..   In  the  next  section^  we  study  this  aspect 
in  greater  detail. 

We  also  observe  that  P^  (and  also  3^^,,   P^;  etc.  if 
we  had  continued  on)  will  vanish  at  A  =  1.   This  merely 
says  that  the  equipotential  line  A  =  1  is  also  a  constant 
speed  line?  thus  the  line  A  =  1  and  the  line  A  =  1  coincide , 
as  can  be  seen  already  from  equation  (9)-> 

In  section  5,  we  shall  give  some  numerical  examples,, 
and  also  give  some  more  physical  interpretations  of  these 
results, 
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A|..   MORE  DETAILED  INVESTIGATION  OF  THE  EQUATIONS  AND  THE 

SINGULARITIES 

The  basic  equations  (11)  -  (14),  while  convenient  for 
calculating  the  power  series  expansions,  do  not  reveal  the 
nature  of  the  singularities.   We  re-write  (11)  and  (12) 
(ignoring  the  harmless  and  uncoupled  equations  (13)  and 
(14)   ),  using   dp  =  c  dp   ,  into 


o  '  ^ 


The  characteristic  directions  are  given  by 


dO  ^  _   (A^-1)^  (M^-1)     An^i! 


h 


^^M   "(a/.1)2  (m2-.a2.1)   ">  q^2  ^24) 

2   2 
so  that  the  axis  becomes  characteristic  at  M  +A  =  1,  and 

(seemingly)  also  at  A  =  1,   The  sonic  transition,  M  =  1, 

offers  no  difficulty,  since  there  the  axis  is  perpendicular 

to  the  characteristic  directions,  and  hence  is  noncharacter- 

istic.   We  handle  each  singularity  separately.   Fig.  2(a) 

shows  the  characteristics  in  the  nozzle  flow. 
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(1)  M^+A^  =  1 

The  characteristics  here  are  tangent  to  the  axiSj,  and 
the  transition  line  is  perpendicular  to  it.  To  illuminate 
the  situation;,  we  consider  the  Tricomi  equation 

xu   +  u   =  0  (25) 

yy    XX 

which  has  characteristics  as  shown  in  fig.  a(b).  The  x-y 
axes  are  the  counterparts  of  the  §-ti  axes  in  our  problem, 
the  origin  corresponds  to  our  singularity,  and  the  y-axis 
is  the  transition  line  (although  m  our  problem  the  trans- 
ition line  is  not  the  T^-axis,  but  is  only  tangent  to  it  at 
the  singularity)  o   Cauchy  data  an  the  x-axis,  i.e.  u(x,0)- 
and  u  (Xi,0)j,  cannot  be  arbitrarily  posed  at  x  =  0„ 

We  restrict  our  attention  only  to  solutions  symmetric 
about  the  x-axis„   Since  the  y-axis  is  non-characteristic j 
we  may  consider  every  analytic  solution  in  the  neighborhood 
of  the  origin  as  the  (unique)  solution  of  a  Cauchy  problem 
with  analytic  data  u(0,y)  and  u^lOj,y)  on  the  y-axis.   Thus, 
if  we  give  arbitrary  analytic  data  (symmetrically)  on  the 
y-axis,  we  shall  be  able  to  determine  the  form  of  every  symmetric 
analytic  solution  in  the  neighborhood  of  the  origin,  and  in 
particular,  its  behavior  on  the  x-axis  at  the  origin.   This 
will,  in  turn,  enable  us  to  pose  the  right  kind  of  Cauchy 
data  en  the  x-axis » 
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Thus  J,  giving 

/    \  2      4 

°    ^      ^  (26) 

u^(0,y)=  b^  +  b-^y  +  h^y  „   .   . 


substituting  into  (25)  a  power  series  solution  satisfying 
(26)  and  computing  coefficients  term  by  termj,  we  find  on 
the  X-axis, 

u(x,0)  =  a^  +  b^x  =  (a^/3)x^  =  (b^/6)x^  +  (2a2/15)x^  =-= 

(27) 

and  J  naturally  J  u  (xpO)  =  0 

The  form  of  (2?)  suggests  that  the  data  u(xjO)  to  be 
posed  on  the  x-axis  should  be  a  power  series  with  the  2nd, 
5th5  8th,  11th,---  powers  missing;,  thus  3 

u(x,0)  =  A^  +  A^x  +  A^x^  +  A^x^  +  A^x^  =--       (28) 

Indeed,  if  we  take  any  convergent  power  series  of  this 
form  as  data  on  the  x-axis,  and  repeat  the  formal  procedure j 
we  shall  find  a  formal  solution  u(x,y).   On  the  y=axis, 
this  formal  solution  u(o5y)  and  its  derivative  u  (0,y) 
indeed  do  converge  (with  a  smaller  radius  of  convergence 
than  u(x,0)  ).   Therefore,  applying  the  Cauchy-Kowalewski 
Theorem,  the  formal  solution  u(xj,y)  is  convergent,  and  is 
an  actual  solution,  in  the  neighborhood  of  the  origin. 
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This  proves  that  Cauchy  data  in  the  form  of  (28)  are  ad- 
missible across  the  origin. 

The  same  procedure  will  now  be  applied  to  our  problem. 
Let  the  origin  of  the  i-'^      plane  correspond  to  the  inter- 
section of  the  first  transition  with  the  axis.   We  give 
arbitrary  analytic  data  on  ^  =  0,  and  determine  the  behav- 
ior of  the  (analytic)  solution  on  ^l  =  O;  this  in  turn  gives 
us  the  form  of  analytic  data  to  be  posed  on  T)  =  o,  as  is 
needed  in  the  problem.   Specifically,  we  want  the  behavior 
of  h(  i    ,    0),  which  is  needed  in  (20)  and  (21). 

As  before,  for  any  function  f(  i ,  "^  )  ,    f^  denotes 
f  (  ^  ,  0).   Corrospondingly  f°  shall  denote  f(0,Ti  ). 
Evidently  f°  =  f(0,0)  =  f^   . 
We  give 


q°  =  q^-  (l-^a^-n^+a^T)^-  .  .  ,  ) 


(29) 


0°  =  b^T^b^Ti^+b^n-^ 


and  seek  for  solutions 


o  /  . 
q  =  q   (  X 


-A,;'  ^  A,^^  ^  A..?,^  +  .  .  .  ) 


(30) 


0=0   (1+B^  +  B    e    +  B.^e"  -^  .  .  .  ) 


Where  A^ ,  k^,    ,  B^,  B^,  are  functions  of  t]  alone, 

and  are  to  be  determined. 

In  what  follovjs,  we  shall  let  the  adiabatic  constant  Y 
be  3/2.   This  simplifies  the  algebraic  manipulations 
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considerably.   The  conclusions  reached,  however,  are  also 
valid  for  other  values  of  y  .   Expanding  all  the  other 
variables,  we  have 


1  -  ^  M°2  [i  .  (a_)2j 


Y-1 


1  -  M°2  A^^   +  ^2^2  +  c  p3  + 


h 


h^ 


Q.      P 


3-  ^  =   1  -  D^^  -  D2?2  _  0  j3 


(31) 


2'    ^   -      '^    -      u  '  ^  -  |l  -  ~^    >   .0 


(1  _  ^)  q  =   (1  .  a_J)  q  = 
A 


P    J 


=  (1  - 


-fe)q°  [l  +E^e  -E2^2  ^g^,3  , 


(1- — o)qn   (1-  -?)qi 


A^I 


(1--^)  q    ( 


,02''^ 


-1 


=   G  [l-E^e  +  (Ego-Es+E;^^)^^  +  (^30-^3  ""  2E1E2 
where:  ~'^1^20"'^1  ^^  "*"•• 


^2  =  IT  \  -  -^  ^^1  ■"  2A2) 


M 


o4 


^3  -  4  "iv"i 


At(At2  +  2A2)  -  M°2(A^  A^  +  A-) 

«o^ 


^^  "  TS"  ^^1   "^  2A2)'^  +  ^  A^(A-^A2+A3)  -  ^  (A2'^+2A^A<^+2A^ ) 


D^  =  (1  -  M°2)  A^ 

D2  =  (1  -  M°2)  A2  -  A^2(;L  -   2 


^°'  +  ^  M°^ 


(32) 
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3 


A^(l    -   M°^)    -    2   A^A2(1    -  -^ 


+  ^  M°^)+A^3(l    - 


j^o2        ^06 


=    0 


\ 


r° 


^2^  ,o2  ^ 


A""-    1 


A.2   h02  m^_    ^  W^.l^l\ 

^    2   "     4    ^    '      20"   ^  2_   ^^2'"^^ 


A°2_    1 


E^   = 


,    ,    M°2    .    A°2    -    1    .    ^   V2^    (, 

^3^  — n — r~  ^  ^  "7^2       ^^ 


M 


o2 


G      = 


A,  ^    A^^   M- 


M 


-    1 
2 


■)  + -^dJ^'^^)  ; 

2(A°2-1) 


^    _  ^10  ^20  '•'!    ,0    ':i_) 

'^  a/    -   1 

Al 

(1    -   ^)q° 
A°^ 


''^'    ""lO^   -    (1    -   m2) 


2(A  2    _  1) 


Substituting  (30)  and  (3I)  into  (11)  and  (12),  and  equating 
the  coefficients  for  each  power  of  ^  ,  we  obtain: 


B,   =   - 


h°(l-M°2) 
1  dG 

h°  e°      ^^ 


dO_ 

dn 


B, 


2h°Gg         d^,         1  2  11 


'^  -'ir^  [ 


d  r^      ^Ov  py^         d0 

in-  ^^1'  ^  -  ^"1  r^ 


o 


\ 


(33) 


(34) 
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°3  =  -^ 


G^(B^e°)  -OE^^irf')  *0(E^,.E2*E^^)§ 


df] 


2xd&: 
i 


We  are  after  the  behavior  of  h(  |  ,0);  hence  we  must 
find  D  Q  and  D^^^.   V/ith  q°  and  0°  given  by  (29),  we  have, 
from  (20), 


30' 


G  =  1  -  a^( r-2 —'^     -  \a,,{ ^5 )   -  a^  ( ^(^_-i_). 


■a~2"  _  1 


M  2   M  2 


4''    2 

^   Aj"^  -  1 


1-A^ 


..  V 


Thus 


=  1  "2<3--r)(7^)«2'n*'+---- 


-  2^I°20 


G  5^  (B,e°)  .  GE,  1^ 


T1=0 


as  anticipated.   Now 


'dO^ 


^'-t^H-.:'"^A-i_:° 


all   other   terms   vanishing  since  E,„   =  0, 


Using  B^   from   (33), 


(35) 


(36) 


'I''30  ~  I  dT)    ^  "F 


-  6h,D,.  =  1:3^  (-tc) 


^  (^1^^ 


Tl=    O 


r]=    0 


^tfc'^^^y 


i=    o 


+    (1-M/)    a,,    [^  B,0 


I    '    ""10    \dr\  "1 


+    IB^O 


Tl=     O 


01  /J_(i.M°2)A, 


dr) 


But 


d|-<=i°' 


=      E 


^=  o  "        ^0   I^Vti=  o 


dG] 


r)=   o 
/dE 


T]=     O 


+         1 


1^.1= 


d_    .    (M°^.Af   -   1  , 


dT]      "1^        ^o2  _  ^ 


/A  (_g.)(MZ)  m°2+a°2-i 


dTi      ,0'^dTi'    ,^02 


h 


(M^^-1)(A°2.i)- 


=      0 


T1=     O 
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ri=o 


since  the  quantity  being  differentiated  is  an  even  function  in  n 
ilarly,  /   2      I         / .2o\        dE 


bim- 


^  '^^1 


=  E 


r\=   o 


dTl^ 


d^G' 


+  2(— ^)(^) 


+  1 


d^E^  \ 


r\=   o 


G  dO 


G      dO^  (tl^itA!^)" 
h°(M°2_i)  d^'    A°2  _  1 

2 


o' 


\h 


o  dn 


JJ_  ^  M°^+A°^-l 


^Tl  =  O 


dTj' 


{M°2_i)(A°2_i) 


n=  o 

-)) 


'n=  o 


( 


M°2+A°2_i 


With  q°  given  by  (29),  and   p°,  c°   in  terms  of  q"  given  by  (19) 
and  (20)  respectively,  we  readily  find  that 


M°2  +  A°^  =  Mj^  +  A-j.^ 


M- 


(M,2a^2  +  (2  +  _L_  )m  2)^  n 


1  "I 


and  that 


M°2+A°2_i 


dT^    (M°2.i)(A°2_i) 


iT=  o 


(M-|.2.l)(Aj2.i) 


-%(M°2^A°2) 
dT        / 


(Mj2_i)(A^2_-^) 


(Aj'^Mj'^+(2  +  -|-)Mj  )2a2  =  0 


Thus,  D_^  is  non-zero,  and  is  moreover  independent  of  D^^  (which  is 


r      n' 
equal  to  I(M°2.i)a       =  G^  4^ 
L       -^  Jt>=  o   Ih    'il 


) ,  since  b-,  and  a^ 


are  both  arbitrary  and  thus  independent.   Thus,  we  see  that  the  data 
to  be  given  on  the  i     axis  must  be  of  the  form 

h(  ^,  0 )  =  h^  +  D  ^  +  D    e  3  + ( 37 ) 
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For  h   of  this  form,  equations  (20)  and  (21)  will  then  give  us  the 

first  and  second  approximations  of  the  flow  streamlines,  without  any 

2     2 
further  difficulty  from  the  singularity  at  M   +  A  =  1  . 
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We  shall  content  ourselves  with  choosing  third  order 
polynomials,  for  h  ,  since  investigation  of  higher  order 
terms  becomes  algebraically  prohibitive.   A  third  order 
polynomial  h  will  give  a  flow  crossing  all  three  trans- 
itions, since  such  a  function  can  first  decrease  and  then 
increase,  thus  giving  convergent-divergent  nozzles.   All 
higher  order  coefficients,  moreover,  can  also  be  determined, 
in  principle,  and  will  be  finite  everywhere.   Unlike  the 
case  of  the  Tricomi  equation,  unfortunately,  we  have  no 
proof  that  if  h  is  a  convergent  series,  the  resulting 
formal  solution  will  also  be  convergent,  and  thus  be  a 
true  solution  in  the  neighborhood  of  the  origin.   This 
gives  yet  another  reason  why  our  solution  remains  a  formal 
solution. 

If  we  had  taken  higher  order  terms  in  the  formal 

expansion  given  in  section   3,  we  would  obtain,  as  the 

next  higher  term  corresponding  to  (21)    ^      ^ 

1       A  2  f     n  ^       „     „^  nj         a/  -  1 

p,,    .    1  _o_^_L.      _      o_  h   2   h    2    (       _   _2_)  Q 


^   ~   ^   M    2^A   ^.1      r      8        o      ''o       ^'    -      2     '   ,^    2^_^   2_^)2 
o        o  \_^  00 

where      6^    is   given  by    (21).      Near   the   singularity, 

— M   "   ^ tends   to  a   finite   limit.      Hence   the   second 

Ao2+Mo2  -  1 

and  third  terms  in  the  bracket  tend  to  zero.   If  h  has 

o 

been  chosen  according  to  (37),  the  first  and  last  terms  in 
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the  bracket  will,  undoubtedly,  have  their  zero  order  terms 
( in  ^  )  cancel,  so  that  the  entire  expression  in  the  bracket 
starts  off  v/ith  first  power  in  I  .   Thus  P  ^  will  be  finite. 
However,  direct  verification  of  this  is  difficult  because 
of  extremely  cumbersome  algebraic  manipulations. 

R.  Seebass   has  obtained  a  special  exact  solution 
across  the  first  transition  (hypercritical  transition,  in 
his  terminology)  for  y  =  2.   The  resulting  h^(x)  from  his 
solution  is  exactly  of  the  form  (37).   Kogan   has  also 
given  a  special  solution  across  this  transition.   His 
h  (x)  has  an  x  -term,  which  seemingly  negates  our  result, 
Kogan' s  solution,  however,  is  approximate,  and  is  valid 

only  to  linear  terms  in  x.   Hence,  no  information  may  be 

12 
extracted  from  it  for  the  behavior  of  higher  order  terms. 

(2)  A  =  1 

Equation  (2^)  seems  to  indicate  that  the  axis  is  again 
characteristic  at  A  =  1.   This,  however,  is  not  the  case. 
Consider  the  line  on  which  A^(^)  =  1.   This  is  an  equipo- 
tential  line.   But  (9)  shows  this  also  coincides  with  the 
line  A  =  1.  this  fact  has  been  observed  earlier  at  the  end 
of  section   3.   Therefore,  in  (2^),  the  characteristic 
direction  is  d^/dn  =  0  at  this  singularity,  which  implies 
that  the  axis  is  perpendicular  to  the  characteristic 
direction,  and  is  thus  non-characteristic.   The  data  given 
on  the  axis,  therefore,  should  not  require  any  special 
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handling  as  needed  at  the  other  singularity  M  +  A"^  =  1  ^ 

However  J,  as  a  result  of  A  =1  (and  A  =  1  simultaneous- 
ly) ^  (22)  and  (23)  indeed  fail  to  satisfy  the  hypotheses 
of  the  standard  Cauchy-Kowalewski  Theorem  at  this  singular- 
ity.  For^  in  applying  the  Cauchy-Kowalewski  Theoreraj  A 
must  be  considered  as  a  function  of  the  independent  vari- 
able E,  5  while  A  is  considered  as  a  function  of  the  depend- 
ent variable  q,  and  the  derivatives  with  respect  to  Y)    must 
be  solvable  from  the  equation  Vi/ithout  recourse  to  the  fact 
that  A(  ^  J,  0)  =  A  ;  i.e^j,  the  equation  must  have  analytic 
coefficients  and  be  non-characteristic  along  the  initial 
line  independently  of  the  Cauchy  data.   Nevertheless j  the 
coefficients  of  the  formal  expansion  are  finite  and  apparent- 
ly reasonable;  this  leads  us  to  believe  that  the  singularity 
is  actually  harmless  in  nature j  and  does  not  affect  the 
validity  of  the  solution.   A  refined  mathematiciil  investig- 
ation of  the  well-posed  nature  of  the  Cauchy  problem  near 
this  singularity  will  not  be  attempted  in  this  paper. 
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5.  NUMERICAL  EXAMPLES. 

We  shall  now  give  two  examples  of  flows  calculated  by  the  method 
presented  in  the  previous  sections.   Since  we  are  interested  in  math- 
ematical properties  rather  than  practical  usefulness  in  this  paper, 
we  take  3rd  order  polynomials  for  h   as  required  in  Section  4  to 
insure  well  behaved  solutions*  the  nozzles  obtained,  however,  are 
rather  long  and  not  too  practical  for  actual  use. 

Nozzle  "A"  is  calculated  with  given 
h^  =  5  -  64  +2  4^ 

and  nozzle  "B"  with 

h^  =  3  -  34  +   4^ 

Once  we  have  h  ,  and  the  fact  that  M^^  +  A^^  =  1  at  4  =  0,  all  the 

O '  0       0 

other  flow  quantities  on  the  axis  are  completely  known.   In  particular, 
we  easily  obtain  for  Y  =  1»5, 

which  permits  the  determination  of  M^( 4  ).   For  4  =  0,  we  then  have 

2      2  ? 

M,  and  A-.  =  1-M-r  .  But  since  A  is  inversely  proportional  to  p,  we 

immediately  have  >  _l_ 

'l^^k   m2  1  Y-1 
?      ?        2 


We  shall  confine  our  interest  to  the  region  -0,5  :<  4  ^1.5, 
which  is  the  length  of  the  nozzle.   At  4  =  1.5,  the  flow  is  super- 
sonic (hyperbolic),  and  the  nozzle  contour  beyond  that  can  be  easily 
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shaped  by  simple  wave  theory.   On  the  other  hand,  at  the  inlet,  we 
do  not  have  a  flow  at  rest,  but  only  a  very  low  speed  flow  with  a 
given  speed  distribution  across  the  inlet  section.   That  this  can 
be  physically  realized  we  shall  assume. 

We  also  observe  that  when  the  flovj  is  at  rest  in  the  reservoir, 
M  =  0,  but  A  ^  0,  since  A'  =  — ^   .   Vie  denote  this  by  A   .   For 
each  nozzle,  we  have  the  following  constants: 

Nozzle     M   inlet    M   exit    M,    A^.     A^ 
o  o  II 

"A"        0.074       2.66    0.123   .987   0.980 
B"        0.131       2.19    0.195   0.962   0.943 


lip  II 


Using  equations  (20)  and  (21),  the  flows  are  computed  and  shown  in 
figs.  3  ai^d.  4  respectively.   The  first  approximation  (hydraulic 
theory)  is  superimposed  on  the  same  figure.   The  close  agreement 
between  the  two  approximations  tends  to  suggest  probable  convergence 
of  the  formal  power  series.   In  fig.  3  (Nozzle  A)  streamlines  ere 
shown  up  to  Tj  =  0.1,  and  in  fig.  4  (Nozzle  B),  up  to   ^  =  0.2   This 
means  that  the  throat-length  ratios  are,  respectively,  10;^  and  20;^, 
a  rather  long  nozzle  in  each  case.   Larger  values  in  t^  result  in 
wide  difference  between  the  second  and  first  approximations,  suggest- 
ing that  higher  order  terms  must  be  included,  or  even,  perhaps,  that 
the  formal  solution  is  no  longer  convergent. 

The  positions  of  the  transition  lines  are  of  interest.   They 
are  also  shown  in  figs.  3  and  4.   Their  locations  can  be  qualitatively 
seen  by  the  simple  considerations  as  follows:   The  line  A  =  1  is, 
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as  mentioned  before j,  an  equipotential  line  as  well  as  a  constant- 
speed  line.,  so  its  location  is  obviously  coincident  with  ^  =  ^^-j.. 
In  the  vicinity  of  the  line  M'"  +  A^  ==  1>  or  i  =   iy,    the  quantity 
Pp  is  negative  0  This  implies  that  in  this  vicinity,  the  speed  de- 
creases away  from  the  axis  on  an  equipotential  line.   A  constant- 
speed  line^  therefore^  must  bend  forward  (i.e.  downstream)  of  an 

2   2 
equipotential  line.   Hence  the  transition  line  M  -^A  -  1  bends 

dowstream  from  the  equipotential  line  ^^  ^   0.   An  identical  argu- 
ment shows  that  the  sonic  transition  M  =  1  bends  backwards  (or  up- 
stream) away  from  the  axis;  this  is  the  same  as  in  ordinary  gas 

nozzles . 

We  conclude  with  a  remark  comparing  cur  results  with  the  flow 
of  an  ordinary  non-conducting  gaso  The  latter  can  be  computed 
either  by  the  procedure  of  Friedrichs  for  plane  flow,  or  simply 
by  setting  the  magnetic  field  to  zero  in  our  results  (setting  a  =  0) « 
We  have  already  observed  that  the  first  approximation  (hydraulic 
theory)  remains  the  same.   For  the  second  approximation  (21),  Pg 
simply  becomes  h  h  "/2.  y,  changes  accordingy.*  while  x^   remains 
the  same.   In  our  examples,  the  flows  will  not  be  widely  different, 
since  the  second  approximation  is  close  to  the  first, 

A  more  detailed  study  of  the  effect  of  the  field  on  the  flow, 
involving  more  numerical  calculations,  as  well  as  the  more  practical 
aspects  of  these  flows  (e.g.  comparison  with  other  field  orienta- 
tions, practical  reallzabllity ,  etCc),  will  be  considered  in  a  sub- 
sequent worko 


-  28  - 


ACKNOWLEDGMENTS 

The  writer  is  indebted  to  Prof,  H,  Grad  for  suggest- 
ing the  importance  of  this  probleiru   He  also  wishes  to 
thank  many  of  his  friends  and  colleagues j,  m  particular ^ 
Profs  o  Ko  Oc  Friedrichs,,  C„  S,  MorawetZ;,  and  W„  R^  Sears, 
for  stimulating  and  fruitful  discussions,,   The  numerical 
computations  were  competently  carried  out  by  Mr.  Y.  Wu. 
This  research  was  initiated  while  the  author  was  at  the 
Institute  of  Mathematical  Sciences  during  the  summer 
of  I96O5  and  was  supported  by  the  Uo  S„  Atomic  Energy 
Commission  under  contract   AT  (30  -  1)  1480  5  it  was 
concluded  during  the  academic  year  I96O  -  61  under  a 
National  Science  Foundation  grant  G  1^1712  0 


29 


FOOTNOTES 


1.  H,  Grad,  'Reve  of  Mod.  Phys.,  32,  830  (I96O). 

2.  K.  0.  Friedrichs,  "Theoretical  Studies  on  Flow  Through  Nozzles 
and  Related  Problems",  Nat,  Defo  Res.  Comra, j  Appl,  Math.  Panel, 
Report  82.1  R  'a9J^J^). 

3.  F.  I.  Frankl,  Izvestiya  Akad.  Nauk  SSSR,  Ser,  Matem. ,  9,    387  (1945; 

h.        The  situation  of  more  equations  than  unknowns_is  frequent_ly 

encountered  in  physical  problems,  e.g.,  curl  q  =  0,  div  q  =  0 
in  classical  hydrodynamics. 

5.  M.  N.  Kogan,  Prik.  Mat.  Mekh. ,  23,  70  (1959). 

6.  E„  L.  Resler.  Jr.,  and  W.  R,  Sears,  Jour,  of  Fluid  Mech., 
5,  257  (1959)0 

7.  If  the  field  and  velocity  are  always  in  one  plane  (two-dimensional 
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velocity  and  field  everywhere. 
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9.  We  could  just  as  well  have  considered  axi-symmetric  flow,  with 
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more  useful  result.   Nevertheless,  the  essential  mathematical 
features  are  no  different,  and  we  shall  do  the  simpler  case, 
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FIG.  2(a).  Characteristic  curves 
In  nozzle  flo.v. 
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FIG.  2(b).   Characteristic  curves 

for  Trlcoml  equation  xu  +  u   =  0. 
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